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The most general solution to the Einstein equations in 4 = 3 + 1 dimensions in the asymp- 
totical limit close to the cosmological singularity under the BKL (Belinski-Khalatnikov-Lifshitz) 
hypothesis, for which space gradients are neglected and time derivatives only are considered, can 
be visualized by the behavior of a billiard ball in a triangular domain on the Upper Poincare Half 
Plane (UPHP). The behavior of the billiard system (named 'big billiard') can be schematized by 
dividing the succesions of trajectories according to Poincare return map on the sides of the billiard 
table, according to the paradigms implemented by the BKL investigation and by the CB-LKSKS 
(Chernoff- Barrow- Lifshitz- Khalatnikov- Sinai- Khanin- Shchur) one. Different maps are obtained, 
according to different symmetry-quotienting mechanisms used to analyze the dynamics according 
to the symmetries of the billiard domain and to the features of the geodesies on the UPHP. In the 
inhomogeneous case, new structures have been uncovered, such that, in this framework, the billiard 
table (named 'small billiard') consists of 1/6 of the previous one . The properties of the billiards 
described in these cases have already been studied in the reduced (according to the Poincare return 
map) phase space. 

The connections between the symmetry-quotienting mechanisms are further investigated on the 
UPHP. The relation between the complete billiard and the small billiard are also further explained 
according to the role of Weyl reflections. 

The quantum properties of the system are sketched as well, and the physical interpretation of 
the wavefunction is further developped. In particular, a physical interpretation for the symmetry- 
quotienting maps is proposed, according to the Fourier decomposition of the energy levels of the 
wavefunction of the universe, as far as the meaning of epochs and eras is concerned from a quantum 
point of view. 

PACS numbers: 98.80.Jk Mathematical and relativistic aspects of cosmology- 05. 45. -a Nonlinear dynamics 
and chaos; 
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I. INTRODUCTION 



The description of the chaotic behavior of a generic cosmological solution was approached in Q, 0, g[, 0], @ 
by locally reducing the corresponding Einstein equations to a system of ordinary differential equations (ODE's), 
and then by finding an approximated solution for the system of ODE's. In particular, this reduction is possible 
under the assumption that, close to the cosmological singularity, spacial gradients can be neglected with respect to 
time derivatives, such that points are spatially decoupled. The resulting behavior is one consisting of a succession 
of (suitably - joined) Kasner solutions. The statistical analysis of the model was achieved in [l| - Q, where the 
one-dimensional BKL map (named after Belinskii, Khalatnikov and Lifshitz) was established. This statistical analysis 
was given a precise geometrical interpretation in , 0] , @ where the statistical properties of the BKL map were 
connected to the features of a billiard on a Lobachesky plane. 

In 03, [12 and [13, the CB-LKSKS (named after Chernoff, Barr ow, Lifshitz, Khalatnikov, Sinai, Khanin and 
Shchur) map is defined by upgrading to BKL map for a system of two variables. The BKL map is obtained from the 
CB-LKSKS map by marginalizing the 'extra' variable in the definition of probabilities associated to the map. The 
definition of all these maps is based on a suitable symmctry-quoticnting mechanism. A link between two-variable 
maps and the discrete map for billiards in the Lobacesky plane was given in [l3j . 

The precise relations between the dynamics of the asymptotic Bianchi XI universe and its symmetry-quotiented 
versions (according to the identification of the role of Kasner exponents, and its symmetry-quotiented version ac- 
cording to the geometrical features of the billiard table) were established in [ijj]- This was possible by means of the 
identification of the geometrical expression of probabilities for billiards (given as the integral of an invariant reduced 
form over the pertinent domain in the reduced phase space) and the statistical probabilities for the statistical maps. 
Full information about cosmological billiards from different viewpoints can be found in [l5j , fl(| , fl7j . 

In this work, the features of the unquotiented billiard dynamics and those of its two symmetry-quotiented ver- 
sions are defined for complex variables and analyzed according to their structure, implemented on the Upper Poincare 
Half Plane (UPHP). In particular, the mechanisms of symmetry-quotienting the roles of the Kasner exponents is 
demonstrated to directly imply the proper boundary identifications of the group domains that generate the billiard 
maps on the UPHP. The maps for the reduced phase space given in fl4| are immediately recovered as the limit on the 
absolute of the UPHP. The reduced-phase-space method [14[ and the UPHP one arc illustrated to be complementary 
for the definition of epochs and eras and their succession. 

The paper is organized as follows. 

In Section UH the features of cosmological billiards are recalled. In particular, the main steps in the definition of 
cosmological billiards are outlined, and much attention is paid to the definitions of the unquotiented big billiard and 
of its two possible symmetry-quotiented mechanisms, i.e. according to the permutation of the Kasner exponents and 
to the symmetry properties of the big-billiard table. 

In Section IIII1 the basic tools in the definition of the groups are introduced, out of which the billiard maps are 
implemented in the UPHP for complex variables. 

In Section ITVl the properties of these three kinds of billiards are analyzed on the UPHP. More in detail, a difference 
is found between quotienting out the role of the two oscillating Kasner exponents during a given BKL era, and 
identifying the 'newly-oscillating' Kasner exponent of the next BKL era with one of the previous BKL era. This 
difference implies a difference in the boundary identification for the fundamental domain where the two maps act. As 
a result, a paradigm to construct the unquotiented billiard map, the BKL map and the CB-LKSKS maps starting 
from the geometrical symmctry-quoticnting of the unquotiented domain is established. 

The quantum regime is analyzed in Section [Vj where the properties of the wavefunction, resulting from the solution 
of the Hamiltonian constraint, are fixed according to the different physical features of the billiard table. More in 
detail, the symmetries of the billiard table and of the dynamics are taken into account, as far as the definition of the 
mathematical features of the wavefunction is concerned. 

The role of (BKL) epochs and eras in the quantum version of the model is analyzed in Section fVTl where the meaning 
of the 'quantum' BKL numbers is also investigated. Because of the features of the quantum version of BKL epochs 
and eras, open issues are underlined, and a comparison with previous results of the investigation of cosmological 
quantum billiards on the Unit Hypcrboloid (UH) is sketched. 
Some brief concluding remarks in Section [VIII end the paper. 



3 



II. BASIC STATEMENTS 



The metric of a general inhomogeneous space-time close to the cosmological singularity [l5l | is suitably parametrized 
by the Iwasawa decomposition of the spacial metric, vanishing shift functions N l and an appropriate choice of the 
lapse function N . Close to the singularity, the limit of the Iwasawa decomposition of the spacial metric reads 

d 

giJ (x°,x l ) =J2e- 20a ^°^N a i {x o ,x i )N a j (x°,x i ). (2.1) 

a=l 

The asymptotical limit close to the singularity is characterized by the f3 a functions only, while the remaining degrees 
of freedom (given by the upper triangular matrices A/" a j(a; , x 1 )) are frozen. The corresponding model is that of a 
masslcss particle moving in the (/3) space with metric 



da 2 = J2 G ab d(3 a df3 b = (df3 a f - 

a,b=l a \ a / 



(2.2) 



in a closed domain delimited by infinite sharp walls. The dynamics is therefore that of a billiard, i.e. characterized 
by (elastic) 'bounces' against the walls and 'free-evolution' (geodesies) trajectories between the billiard walls (given 
by the asymptotic limit of the corresponding potential). 

The free evolution between two bounces is given by the Kasner solution, the 'curvature' walls describe the asymptotic 
limit of the homogeneous Bianchi IX model, while the 'symmetry' walls have to be considered for inhomogeneous 
models (details will be given in the next paragraphs). 

The Kasner metric is obtained for a diagonal metric with gij = a(t)dx + b(t)dy + c(t)dz, where a, &, c are the Kasner 
scale factors, such that a(t) = e~" x (and analogous identifications for the other variables). It admits the solution 
pa _ v a T _|_ w } iere t} ie time variable r is defined via dr = —dt/^Jg. The Kasner velocities v a define the Kasner 
parameters p a as p a = v a /J2 a v<1 ■ wmcn obey the constraints J2 a P a ~ Ci-P") 1 = 1- 

When the potential is taken into account, the complete behavior close to the singularity is given by a piecewise 
succession of Kasner solutions 'glued' together at the bounces against the billiard walls, where the role of the Kasner 
exponents is permuted, according to the BKL map as a function of the variable «bkl- 

Decomposing the variables f3 a as f3 a = pj a , such that 7 Q 7 a = —1, solving the Hamiltonian constraint for lap is 
equivalent to project the dynamics on the UH j a "f a = — 1. 

By means of appropriate geometrical transformations (for which no informations about the dynamics is lost), it 
is possible to describe the dynamics on the UPHP [l8| [l^] [20l |. In the UPHP, with coordinates (u,v), geodesies 
are (generalized) circles (with radius r) orthogonal to the u axis (i.e. centered at (u = uq,v = 0)). They can 
be parametrized with their oriented endpoints on the u axis, say u + and u~ , such that r =\ u + — u~ \ /2 and 
i*o = {u + + u~)/2. Although the complete phase space for billiards on the UPHP is in principle 4-dimensional, the 
existence of a reduced 2-dimcnsional form w(w~,it + ), 

du + A du~ . . 

c^ = — -2, (2.3) 

(w + — u Y 

invariant under the billiard dynamics in the reduced 2-dimensional phase-space («~,w + ) explains the choice of the 



paramctrization of geodesies according to their oriented endpoints as completely defining the billiard dynamics [14 j . 
Furthermore, the geometrical definition of probabilities, given by the integration of the form (|2.3p over the pertinent 
subdomain of the (u~, u + ) reduced phas e sp ace, exactly coincides with the definition of probabilities for the statistical 
maps of the billiard, given in [l(|, (Tl|, |l2j |. 

After the projection on the UH, the expressions of the gravitational walls in term of the /3 variables are equivalent to 
those in term of the 7 variables. The Kasner parameters p a can be expressed as a function of the 7 variables, and, 
on their turn, in terms of the Kasner scale factors a, b, c. The names of the three curvature walls in (|2.4[) below 
defining the three sides of the big billiard are due to the fact that the expression of the 7 variables as functions of 
the coordinates u and v of the UPHP reduce to that of the Kasner parameters p as a function of the ubkl variable 
at u = u + and v = 0. The role of the reflections on the curvature walls is then interpreted as the permutation of the 
Kasner parameters p via the BKL map on the variable ubkl (rigorously) identified with u + . The role of the Kasner 
parameters will then be discussed according to all these relations. 



4 



a. The big billiard In the UPHP, the big billiard consists of the billiard table delimited by the curvature walls a, 
b, c, and are described by reflections laws A, B, C given by 

a: u = 0, Au = —u (2.4a) 
b: u=-l, Bu=-u-2 (2.4b) 

c: ( u+ i) 2 + , 2 = I, Cu = -^ n , (2.4c) 

where the transformations A, B and C act diagonally on u (i.e. independently on u + and it - ). The big billiard 
tabic is sketched in Figure [TJ 

The billiard dynamics can be therefore parametrized according to: i) the 'free-flight' evolution between to consecutive 
reflections on the billiard walls, and ii) the three reflections A, B, C against the three billiard walls a, b and c. 
The (appropriate composition of) the transformations A, B, C define the big-billiard map T, and the reduced form ui 
in (|2.3p is invariant under the map T . The billiard dynamics is therefore characterized by the bounces on the billiard 
walls, and the free-flight evolution is neglected. As a result, a Poincare return map on the billiard walls is obtained. 
The dynamics can be encoded as a sequence of epochs and eras. An epoch is defined as a trajectory (free-flight 
evolution) between any two walls, while an era is defined as a collection on epochs taking place in the same corner. 
Accordingly, 6 kinds of (oriented) epochs are distinguished, according to the walls connected by the corresponding 
trajectory. In more detail, these are ab, ba, ac, ca, be, cb. Correspondingly, 6 kinds of eras are obtained, and are 
named after the first epoch starting the era. 

The reduced phase-space (u~u + ) can therefore be divided into 6 different regions, where the 6 kinds of epochs take 



place. These 'boxes' are named B xy after the corresponding xy epochs; their details arc given in Table II of [14j . and 
they are sketched in Figure [3J where the B xy boxes arc delimited by solid black lines. The succession of epochs is 
described in this plane as a succession of 'jumps' in a epoch-hopscotch game. Probabilities for sequences of epochs 
and eras are defined by the integration of the reduced form to in (|2 . 3[) over the pertinent regions of the (u~u + ) plane. 
During a Kasner xy era, the two Kasner exponents x and y oscillate (thus defining epochs), while the third one, 
say z, undergoes a monotonic evolution. In the next Kasner era, the z exponent takes the place of one of the two 
previously-oscillating x or y (according to whether the previous era contains an odd or an even number of epochs) . 

In addition, the compositions of the big-billiard map T allow one to describe the billiard dynamics of epochs, 
for which only the initial values Up determine both the type and the length (expressed as the number of epochs), and, 
according to them, also the following sequence of eras. The corresponding starting regions for u^, in the (u~u + ) plane 
can be defined, and the dynamics is described as a succession of 'jumps' in an era-hopscotch game. The definitions of 
the F xy regions are given in Table III of [l4| , and they are sketched in Figure [3l where the F xy boxes are delimited 
by (dotted) blue lines. 

b. The Kasner- quotiented big billiard The Einstein equations from which the billiard dynamics is derived are 
invariant under the group of permutations between the three functions a, b, and c. The corresponding billiard 
dynamics can be therefore quotiented by use of this group of permutations by defining BKL epochs and BKL eras. 
During a single Kasner era, the role of the oscillating x and y exponents is identified in a BKL epoch, while that of 
the monotonically-evolving z is kept fixed. The transition to the following BKL era is obtained by identifying the role 
of the z exponent with either x or y, and keeping the remaining one (either y or x, respectively) as fixed. 
This is equivalent to mapping all the starting boxes B xy (with the corresponding F xy ) onto a (preferred) one, say ba 
(for technical purposes). This is obtained by means of the Kasner transformations fc;. The explicit expressions of the 
Kasner transformations fc^'s are listed in [1J], and are obtained from Table U by imposing v = 0; their geometrical 
interpretation in permuting the Kasner exponents, as well as their specific action with respect to the (preferred) 
epoch-type ba, are given in Table HI 

The corresponding quotiented map Tbkl is a succession of compositions of the Kasner transformations ki and the 
big-billiard map T. The usual BKL map for the u = = mbkl variable, Tbkl 

u -> u - I -> u - 2 -> ... -> u — n - I -> u' = m _ 1 j m j — 1, (2.5) 

is then obtained. In (|2.5[) , [u] denotes the integer part of u, and [it + ] = n — 1, n being the number of epochs in the 
u er4E such that < [up\ < n — 1. In (|2.5p . it is possible to identify two different features: i) the oscillatory regime 



1 The BKL u> o definition for the «bkl variable, defined in Eq. (5.4) of [5J and used in [T3 . is used here; as it was discussed in |14|I . this 
definition is more suitable to naturally describe the features of the billiard. Of course, also the BKL u> i can be used, after a proper 
relabelling (rcscaling) of the variables. 
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between two Kasner exponents (say a and b) in the part u — > u — 1 — > with the other Kasner exponent (say c) 
evolving monotonically, i.e. the 'BKL epoch map'; and ii) the transition to the next BKL era, when the c exponent 
takes the place of one of the two previously-oscillating exponents (whose role is quoticntcd in the BKL map), i.e. the 
'BKL era-transition map', corresponding to the v! = — 1 element. 

During a BKL era, the BKL epochs correspond to geodesies with the same radius r, and cantered at v = with uq 
scaling according to (|2.5|) . The next BKL era consists of geodesies of different radius r' and different center u' , v = 0, 
r' and u' Q being defined by the BKL era-transition map. 



It is possible to extract from the BKL map (|2.5[) information about the first epoch of each era only, and ne- 
glect the succession of epochs in each era, considering the boxes F xy only. For the particular choice of the 6a-epoch 
type, the era-hopscotch dynamics is the Chernoff- Barrow- Lifshitz- Khalatnikov- Sinai- Khanin- Shchur map T for 
the big billiard, or CB-LKSKS map in short. This map is defined as 

Tu f = + ± +1 - 1- (2-6) 
u F - [u F \ 

In this case, the properties of the dynamics are encoded in the succession of eras, defined by their first epoch only. In 
particular, the CB-LKSKS map is a succession of BKL era-transition maps. Correspondingly, a Poincare return map 
on the a wall for the first epoch of each quotiented ba era is obtained. 

Both the BKL map Tbkl an d the CB-LKSKS map are based on the identification of the role of different Kasner 
exponents, and therefore destroy the typical 'bouncing' description of billiards, based on reflections on the billiard 
sides. As a result, a BKL era ba is described by a succession of translated ba epochs, where no bounces are considered. 
One of the key-features of these maps is the possibility to extract the exact number of epochs contained in the 
succession of eras from the continued- fraction decomposition of the variable u + ; this feature has allowed for a rigorous 
statistical analysis of the BKL dynamics. In particular, this expression has a very simple functional expression on 
the case of ba epochs and eras. Nevertheless, it is possible in principle to define a BKL map and a CB-LKSKS map 
for each one of the xy kinds of epochs and eras; their explicit form would be more complicated than (|2.5p . and the 
extraction of the number of epochs in each era less immediate. 



c. The small billiard The small billiard is obtained in the more general non-homogeneous case. 
The small-billiard table is delimited by the three walls G, B, R, which stand for Green, Blue and -Red, respectively, 
and whose corresponding reflection laws define the small-billiard t map 

G : u = 0, Gu = -u, (2.7a) 
B : u = -\, Bu = -u-l (2.7b) 
R: u 2 +v 2 = l, Ru=±. (2.7c) 

In particular, The G wall consists of part of the curvature wall a in (|2.4a|) , while the B wall and the R wall are (parts 
of) symmetry walls. The 'color' names G, B and R are inherited from [14j]. The small billiard table is sketched in 
Figure [H where a possible use of the 'color' names is depicted. 

Epochs can be defined also in this case as trajectories joining any two walls, but, according to the geometry of the 
small-billiard table, small-billiard eras are defined as a sequence of epochs between two reflections on the R side. The 
regions of the reduced phase-space (it~, u + ) where the dynamics of the small billiard takes place are described are 



BG : u~ < -I, u + > ii„; -I < u~ < -1/2, u + > up; (2.8a) 

BR : u~ < -1, up < u + < u a ; (2.8b) 

RG : -1 < u~ < -1/2, u a < u + < up; -1/2 < u" < 0, u + > u a ; (2.8c) 

RB : -1/2 < u~ < 0, u + <up; < u~ < 1, u a < u + < up; (2.8d) 

GR :u~>l, u a <u + < up; (2.8e) 

GB : < u~ < 1, u + < u a ; u~ > 1, u + < up, (2.8f) 



where the two functions u a and up are defined as u a = —^p- and up = — 2 U -+i ' respectively. The regions of the 
(u~ , u + ) allowed for the dynamics of the small billiard are classified according to reflection laws on the boundaries of 
the small billiard. In particular, the (u~ , u + ) plane is divided into 6 subrcgions allowed for the oriented endpoints of 
trajectories joining the three different walls G, B and i?, and they are sketched in Figure where the subrcgions are 
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delimited by (dashed) purple lines. The small-billiard is then efficient in reproducing the dynamics of the unquotientcd 
big billiard in the following sense. 

The symmetry walls that define the six 'small billiards' of which the big billiard can be constructed of can be con- 
sidered as 'symmetry lines' of the big billiard itself, as sketched in Figured In particular, the epochs of each era of 
the big billiard always cross the B line, while the crossing of the R lines cannot be defined a priori, as one can verify 
directly. This behavior is in one-to-one correspondence with the dynamics of the small billiard, which shows when one 
of the the two oscillating factors vanishes (the side G of the small billiard is hit, and one of the sides of the big billiard 
is hit), when they equal each other (the side B of the small billiard is hit, and a B line is crossed by a trajectory of 
the big billiard), and when the non-oscillating Kasner exponent crosses one of the other two when changing its slope 
(the R side of the small billiard is hit, and the trajectory crosses a R line in the big billiard). Consequently, each 
era of a given type (say ba) of the unquotiented big billiard can be further (sub)-classified according to whether the 
first epoch and the last epoch of the era cross the R line or not. The sub-classification then consists of five different 
i?-crossing properties for each era. 

The comparison of the dynamics of the small billiard as far as maps are concerned can be obtained by considering the 
Poincare return map on the G wall. The region of the {u~u + ) plane where the dynamics take place shows that small- 
billiard eras contain both ba and ca big-billiard epochs in the RG and the BG regions. In particular, the small billiard 
dynamics reconstructs the complete i?-crossing repertoire for an era, but selects some of them from the B^ a box, and 
others from the B ca . As a result, the small billiard fails in reproducing the BKL map because BKL eras are defined 
according to the Poincare return map on the billiard walls, while small-billiard eras according to the hit on the R wall. 

Nevertheless, a CB-LKSKS map will be constructed for the small billiard, and the meaning of a n-epoch BKL 
era can be recovered. 



To better characterize the features of these billiards, their properties on the UPHP will be investigated. 
The main features of some groups (whose properties will be used in the following) are here briefly revised. For the sake 
of a compact description, only the properties which are relevant in this analysis will be treated, while a systematic 
classification can be found in poj . 

A thorough investigation and a rigorous explanation of group construction of the fundamental domain on the UPHP 
for cosmological billiards is given in [34| . where, also in higer dimensions, the shape of the fundamental domain is 
provided on the higher dimensional UPHP, and an elegant formula for the volume is stated. 

d. The SL(2,Z) group It consists of the two transformations 



III. DOMAINS 



T(z) = z + 1, 



(3.1a) 
(3.1b) 



Its fundamental domain is one delimited by the sides Ax, A 2 , A3, defined by 



Ax 
A 2 
A, 




u = 5 , 
u 2 +v 2 = 1. 



(3.2a) 
(3.2b) 
(3.2c) 



The sides are identified as 



T 



S 



Ax A 2 , 
A 3 -> A3, 



(3.3a) 
(3.3b) 



such that 



Tz = z + l, 



(3.4a) 




(3.4b) 
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e. The SL(2,C) group Its fundamental domain consists of half that of (|3 . 2[) . and coincides with the small-billiard 
table (sketched in Figure [5J black thick solid line) . It is generated by the transformations against the sides 

R 1 (z) = -z, (3.5a) 
R 2 {z) = -z+l, (3.5b) 
R 3 (z) = \ (3.5c) 

, which give the complete z version of the small-billiard map t. The reflection laws in the u~u + plane for the small 
billiard table as a function of the variable u are therefore given by (|3.5p evaluated on the v = axes, where the oriented 
endpoints of the geodesies are defined. Because of the domain of this group is obtained by a desymmctrization of that 
of SL(2, Z), it is not possible to identify sides in any manner. 

It is straightforward to verify that the tranformations (|3.1[) are obtained by a suitable composition of (|3.5[) . i.e. 

T = R 2 Ri, T" 1 = R 1 R 2 , (3.6a) 
S = RxR 3 , S- 1 = S, (3.6b) 

while the converse is not true, it is not possible to obtain (|3.5[) from (|3.1|) . In particular, p.l[) do not include complex 
conjugation, and can be interpreted as a composition of two (an even number of) reflections. 



IV. BILLIARD MAPS ON THE UPHP 



Given an epoch parametrized by the geodesies 

/ n2 , 2 2 U + +U~ U + -U~ 

{u-uo) +v =r , u = , r = , (4.1) 

transformations can be defined for the billiard trajectory parametrized by (|4.1[) . The xy kinds of epochs and eras 
are still defined by the values of u~ and u + , as explained in Section lU defined in Table II and Table III of and 
sketched in Figure [3] 

It is more convenient, nevertheless, to examine the maps on the points z — u + iv (which belong to a given geodesies 
(|4.ip ) that imply the transformation laws for the geodesies (|4.1j) . rather than the transformation laws for the geodesies 
themselves. In fact, the reduction of these maps on the u axis will express the transformation law for the variables 
u + and u~ in (|4.1[) . 

In this respect, it is useful to recall that, in fl4| . the choice to coordinatize the reduced phase space by u~ and u + , 
rather than by u and r, is due to the possibility, for the (u - , u + ) choice, to obtain subrcgions delimited by straight 
lines ('boxes') for the definition of B xy and F xy . 

The interest in this analysis is due to the fact that the complete maps for the z variable in the UPHP automatically 
take into account the way the roles of the Kasner exponents are quotiented. As a result, the proper understanding of 
the symmetry quotienting of the Kasner exponents allows one to impose the correct identifications of the boundaries 
of the billiard tables. These identifications then define the (sub)group structure of the maps. 

In particular, the maps that will be considered are the implementation of the maps T (|2.4[) . Tbkl (|2.5p . T (|2.6p and 
t (|2.7|) for the complete z variable in the UPHP. 

/. The big-billiard map The big billiard is defined on the domain defined in (|2.4a[) and sketched in Figure [TJ No 
identification of the sides is required, as the big billiard map T allows one to follow the dynamics throughout the 
corners, and to keep track of the orientation. 

The big-billiard map T for the variable z = u + iv can be obtained from (|2.4[) , and expressed as a suitable composition 
of the elements of (|3.5[) , i.e. 

z, A = Ri, (4.2a) 
f-2, B = R 1 {R 2 Ri)(R2Ri) (4.2b) 
55^, C = (R 1 R 3 )(R 2 R 1 )(R 2 R 1 )R 3 , (4.2c) 

where, for the sake of compact notation, the original names given in (|2.4p have been kept; anyhow, no confusion 
should arise, as these transformations will act here and in the following on the z variable. 

As a result, one sees that the transformations (|4.2p consist of an odd number of reflections, under which the dynamics 
has to be invariant. 



A : 


z ->• - 


B : 


z -> - 


C : 


z ->• - 
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ca — > ba 


Ki 


z — >• 1/z 


i?3 


ac — > ba 


A" 2 


z ->• -(1 + z)/-z 


(RiR 2 )(RiR 3 ) 


be — > ba 


K 3 


2 -> -Z/(Z + 1) 


(RiR 3 )(R2Ri)R 3 


cb — >• ba 


Ki 


z->--l/(z + l) 


(R 3 Ri)(R 2 Rx) 


ab — > ba 


K B 


2 — >• — Z — 1 




ba — > &a 


K 


Z — i" 2 


i 



TABLE I. The Kasner transformations. 

Right Panel. The action of the Kasner transformations. They act on the set of the Kasner exponents a, b, c by suitable 
permuting them. For the preferred order a, b, c, the Kasner transformations have the effect of mapping xy epochs in ba epochs. 
Left Panel. The Kasner tranformations Ki acting on the variable z = u + iv of the UPHP. They can be expressed as a suitable 
composition of the reflections (|3.5[) . Evaluated at v = 0, the Kasner transformations KiS reduce to the fci's of [3 for the 
variable u and for the reduced phase space variables u . 



g. The Kasner- quotiented big billiard maps The S^-quoticnting of the big-billiard is obtained by defining the 
BKL maps and the CB-LKSKS map on the (it, v) plane for the appropriate z — BKL map. The choice of identification 
of boundaries of the corresponding domain will play a crucial role in the analysis of the maps. 

To define the symmetry-quoticnting maps, one first needs the z version of the Kasner transformations, say Ki. 
The Kasner tranformations Ki for the variable z = u + iv are given in Table U where they are also expressed as a 
suitable composition of the elements of (|3.5[) . For v = 0, the fc, Kasner transformations for the u variable alone of 
Table I and Table VI in [l4[ are recovered. The transformations fc, also act diagonally on the variables (u + 7 u~) of 
the reduced phase space. The determinant of the fej's is equal to 1 for the epoch types oriented in the clock- wise 
sense (i.e. ba, ac and cb, 'parallel' to ba), and to —1 for the epoch types oriented in the counter-clock- wise sense (i.e. 
ab, be and cb, 'antiparallel' to ba). Similarly, the K^s for the epoch types 'parallel' to ba contain an even number of 
reflections (and do not imply complex conjugation), while those for the epoch types 'antiparallel' to ba contain an 
odd number of reflections. This way, as given in Table HI the Kasner transformations K$, Ki and K4 are rotations, 
while the Kasner transformations K\, K% and K$ are reflections. 

The BKL map and the CB-LKSKS map for the variable zbkl is obtained by a suitable composition of the big-billiard 
map (|4.2p and the Kasner transformations Ki given Table Q] for the (preferred &a-type epochs) . 



h. The BKL epoch map on the UPHP The definition of the BKL epoch map can be constructed as follows, for 
the variable Zbkl- For the case of a 3-epoch ba era (which consists of the epochs z\ a , z\ h and z^ a ) , the following 
compositions hold: 

Zba, 4kL = K o4a ( 4 -3a) 

z lb = A (4a) = ~ Z L 4kl = K b A{zl a ) = z\ a - 1 = 4 KL - 1 (4.3b) 

4a = BA(zl) = z\ a - 2, 4 KL = K (z 3 ba ) = z\ a - 2 = 4 KL - 2 (4.3c) 

This way, the BKL epoch map is generated by the (composition of) translations T~ m for suitable values of m, 
< m < n — 1. For the variable zbkl, the BKL epoch map then reads 

ZBKL -> T-^BKL - 1 -+ T- 2 Z BKL -> ... (4.4) 

and, given zbkl = «bkl + m'bkl, one easily obtains 

u — > u — 1, v — > v. (4-5) 



As already specified, the BKL epoch map identifies the the a and b Kasner coefficients and keeps that of the c 
coefficient as fixed. According to the roles of the billiard walls in defining the oscillatory behavior of the a and b 
coefficients by the corresponding reflections A and B in (|4.2p while the c side of the billiard has not been reached yet, 
it is natural to require boundary identifications for the big-billiard table as 



b — > a, z — > z + 1, 



(4.6a) 
(4.6b) 
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As the BKL map consists in describing the dynamics by mapping all the epochs in ba epochs, it is fully consistent 
to remark that (|4.3p consists of an even number of reflections, that compose as translations. The 'bouncing' behavior 
typical of billiards is then completely eliminated. 

It is interesting to remark that the BKL epoch map for the z variable on the UPHP is defined by the bound- 
ary identifications (|4.6[) . which consistently describe the symmetry-quotienting of the roles of the a and b Kasner 
exponents, while the c Kasner exponent is left unchanged, during a given BKL epoch. In the reduced phase space 
(u~ , contrastingly, the BKL epoch map for the variable is defined by the behavior u — > u — 1 — > ... only, for 
which it is in principle difficult to trace back the symmetry-quotienting mechanism. 

Furthermore, the domain with the boundary identifications (|4.6[) takes into account both the quotienting of the roles 
of a and b and the identification of c with itself during the same BKL epoch, in the UPHP for the variable z. On 
the other hand, the BKL epoch map in the {u~ ,u + ) reduced phase space takes into account only for the symmetry 
quotienting of the roles of a and b (while the side c of the billiard has not been reached yet), because this description 
involves only the BKL u variable and not the complete dynamics. 



i. The BKL era-transition map and the CB-LKSKS map on the UPHP The next step is to define the BKL 
era-transition map in the UPHP, and then to express the CB-LKSKS as a sequence of BKL era-transition maps for 
the variable z. 

As in the definition of the BKL epoch map, for the case of a 3-epoch ba era (which consists of the epochs z\ a , z 2 b and 
zl a ) followed by an ac epoch z^ c , the succession (|4.3[) is continued by 

z\ c = ABA(zl) = -z l ba + 2, 4' KL = K 2 ABAz\ a . (4.7) 

(The case of a ba era containing an even number of epochs can be treated analogously and admits the very same 
composition). Eq. (|4.7j) describes the transition to the new era z' , where the role of c exponent is taken by either a 
or b (whose role has already been quotiented). It is straightforward to remark that the BKL era-transition map (|4.7[) 
contains an odd number of reflections, i.e. it implies complex conjugation. For the sake of compactness, this features 
will be discussed after the definition of the CB-LKSKS map. 



As a result, the complete BKL (epoch and era-transition) map for the zbkl variable reads 

z ->• T~ l z -> ... -> T- n+1 z ->z'= - — I = T^SRxT'^z (4.8) 

z — n + 1 

(For the sake of easy notation, the symbols T and S are used in (|4.8|) as a 'shorthand' for their expression (|3.6p as 
functions of (|3.5[) . Furthermore, the subscript bkl has been omitted to keep the notation compact; it will be restored 
when needed to avoid confusion). 

Given zbkl = wbkl + i^bkl, one easily obtains the era-transition map 

(u — n)(u — n + 1 + v 2 ) v , 
u — >■ — . n — z—L, v^- -= ? . 4.9 

(u - n + l) 2 + v 2 (u-n + l) 2 +v 2 y ' 



Neglecting the succession of epochs in each era, the following expression for the CB-LKSKS map T is found 

Tz = I = T^SRiT-^z, (4.10) 

z — n + 1 

which is interpreted as a succession of BKL era-transition maps for the variable z. As analyzed in the above, it is then 
straightforward to recognize that the CB-LKSKS map for the variable z will contain only an odd number of reflections. 

Although the BKL map and the CB-LKSKS do not take into account the bounces on the billiard walls (that 
characterize the unquotiented dynamics), and is, on the contrary, based on a description involving always the same 
(ba) kind of epochs and eras, it is not possible to express the BKL map as containing only an even number of 
reflections at each step. 

In other words, the BKL transition map and the CB-LKSKS map define the new radius r' and the new u value for 
the center u' , v = by means of reflections, although these maps always generate epochs and eras of the 6a-type. 
On the one hand, one sees that the part of the map containing a reflection is the era-transition map, where the Kasner 
exponent c has to be taken into account. The fact that the role of the c Kasner exponent in the BKL map is different 
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form those of a and b could look unexpected, if one considers that the bounces on the three sides a, b, and c of the 
big billiard are described by the laws (|4.2p , which contain all an odd number of reflections- as all bounces are just 
defined as reflections against the billiard walls. 

On the other hand, the reasons for this behavior have to be looked for in the difference between the boundary 
identification (|4.6p and the way the BKL era-transition map (|4.7[) quotients the roles of the Kasner exponents. 
In fact, as already analyzed, the BKL epoch map identifies the exponents a and b, while the c Kasner exponent 
undergoes a monotonic evolution. 

On the contrary, the BKL era-transition map and the CB-LKSKS identify the c Kasner exponent with either a or & 
(already quoticntcd in the epoch map). As a result, these maps have to be defined on a domain with the following 
boundary identification 

b— > a, z^rz + l (4.11a) 

c^a, z^-l + ^— (4.11b) 
z — 1 



which are inequivalcnt to the boundary identification (|4.6I) for the BKL epoch map. From a geometrical point of 
view, the transformation (|4.11b[) is the inversion of the side c with respect to the circle (u — l) 2 + v 2 = 1 (which acts 
a a symmetry line of the big billiard, on which subdominant symmetry walls are defined). 

Given all this, it can be useful to note that, although the BKL era-transition map and the CB-LKSKS map contain 
by construction on odd number of reflections, the intermediate step of the BKL epoch map is a necessary ingredient 
in the definition of the quotiented era maps. Furthermore, it is then possible to identify in the complete BKL map 
the two different parts, i.e. the epoch map containing no reflections, and the era-transition map containing reflections. 

A remark is now in order. As explained in the above, any xy epoch (or era) can be mapped as one of the ba 
type by the suitable Kasner transformation Ki, and the corresponding zbkl will be defined as zbkl = KiZ xy . 
Although the suitable Ki may contain an odd or an even number of reflections (as illustrated in the above, and 
in Table UJ, the considerations developped until here hold for the variable zbkl itself, no matter if in its derivation 
zbkl = KiZ xy an odd number of reflections can be present through Ki. In fact, in that case, the only relevant physical 
information would be that the xy type of epoch is 'antiparallel' to the ba type, without affecting the behavior of the 
map for the BKL variable zbkl- In fact, the choice of the ba type of epochs as a preferred one is merely do to the 
corresponding 'treatable' expressions (|4.8[) and (|2.6[) . rather than to any intrinsic physical meaning. 

j. More about the small billiard on the UPHP As remarked above, and as explained in [l4j |. there is a mismatch 
in the definition of a BKL or CK-LKSKS map for the small billiard in the reduced phase space. 



Anyhow, as remarked in [14j , it is possible to reconstruct the complete map from the reduced phase space by consid- 



ering all the Kasner transformations, given in Table |H and then by considering the following. 

It is important to recall that, in the u + u~ description of the dynamics, trajectories (i.e. epochs) arc identified with 
the oriented endpoints (on the v — line) of the (extended) geodesies, where these trajectories lie. 
The dynamics of each era is described by the behavior of the scale factors. In particular, for a generic era xy in the 
B xy region of the epoch hopscotch court, the two scale factors x and y oscillate between zero and a maximum value 
at each epoch, while the z scale factor (corresponding to the side of the billiard which is not hit during the era) has 
a monotonic evolution. During the xy oscillation, the zero of a scale factor corresponds to the moment when the 
billiard ball collides on a gravitational wall, say x, such that the other scale factor, say y reaches its maximum value, 
and vice versa. During each epoch, by construction, the two oscillating scale factors x and y cross once, when their 
values are equal. In the meanwhile, the non-oscillating scale factor z undergoes a monotonic evolution, and may 
(but not necessarily has to) cross the other two. The beginning of an era is in fact defined as the moment when the 
non-oscillating scale factor z changes its slope, but it can cross the other two scale factors in several ways, as it is 
straightforward to verify. 

As a result, the set of epochs belonging to each eras can be further classified according to these properties. 
The symmetry walls, that enclose the small billiard, can be interpreted as 'symmetry lines' for the big billiard whose 
meaning is suggested by their geometrical construction: each of these three lines represents the points of the billiard 
table, where one of the scale factors equals another. As a consequence, the crossing of a symmetry line by an epoch 
singles out the point of that epoch, where the two scale factors are equal. It is easily verified that an epoch can cross 
one, two or three symmetry lines. 

Each B xy region of the epoch hopscotch court can be therefore further divided into subregions S™, according to 
the behavior of the scale factors, which corresponds to the ordered crossings of the symmetry lines. This regions are 
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S 11 

^xy 


x = y 


o2 

^xy 


y = z, x = y 


q2> 

^xy 


x = y, x = z 


o3 

^xy 


y = z, x = y, x = z 


o3' 
^xy 


x = z, x = y, y = z 



TABLE II. Dynamical subregions of the epoch hopscotch court. 





U~ < — $, U + > U a 

— $ < M~ < —1, U + > U-y 


c2 


— $ < U~ < —1, lt Q < U + < U y 




u~ < -2, < u+ < u a 
—2<u~< — <£>, it 7 < u + < u a 


c3 
>-> ba 


— 2<u~< — $, u 7 < u + < 

— $ < u~ < —1, lt a < u + < Up 


3' 
Sba 


-2 < u~ < -1, < u + < up 



TABLE III. The crossing subregions for the Bf, a portion of the epoch hopscotch court. This Table specifies the details of 
Figure ([5} 



given in Table [TT1 where the different kinds of 'crossings' are explicitly described according to the order they happen, 
and referred to by the index m . 

It is crucial to remark that the five S m regions do not coincide with the partition of B xy into F xy and L xy . 
Let's specify this analysis to one particular region of the hopscotch court, Bf, a , for which the details are listed in 
Table [TTT1 the functions u a and up are defined in (|2.8j) . and u 7 : u + = — "-l^ i the 'Golden Ratio' is approximated 
as $ ~ 1, 618. The two functions u a and up cross at the point u~ = — $ and u + — </>, where the 'Small Golden Ratio' 
<f> is defined as <p + 1 = $. The corresponding subregions for the other kinds of eras are obtained from those given in 
Table Hill by applying the corresponding inverse Kasner transformation of Table H i.e. 5™ = K^ 1 S^. The comparison 
of the dynamics of the small billiard with that of the big billiard is consistent only in those regions of the u + ,u~ 
plane, where it is possible to identify trajectories. In fact, for the big billiard, the u + u~ plane describes collisions on 
gravitational walls. As a consequence, it is consistent to consider only the regions of the u + u~ plane of the small 
billiard, where collisions on gravitational walls are described, i.e. the BG and the RG regions. It is straightforward 
to verify that the BG and the RG regions of the small billiard restricted phase space exactly correspond to the 
regions Sl a , S^ a , S^ a , S^ a , S^ a depicted in Figure [5] and specified in Table Hill By considering the proper Kasner 
transformations, it is therefore possible to recover exactly one B xy region of the epoch hopscotch court, i.e. u — > l/u 
for — 1 < u~ < 0, which allows one to map S^ a , S^ a , S^ a into S% a , S% a , S% a . 

Then, the CB-LKSKS map for the small billiard, tcB-LKSKS, is defined by 

t h2 z = T _1 SiiiT _n+1 «, for(u + ,u _ ) e S£ a and(u+, u _ ) e S% a , (4.12a) 
^',3,3' z = T^SRiT-^Raz, for(u + , u -) e S&,(u+ ir) e S£ a , and(u+ u~) e (4.12b) 

The transformation R3 in the second tcB-LKSKS map is given by the K\ Kasner transformation needed to recon- 
struct the Bb a region in the restricted phase space. 

It is immediate to verify that the two maps contain a different number of reflections. 



As a result, the meaning of the tcB-LKSKS map can be inferred from its action on the description of the dy- 
namics. 

The transformation R3 in (j4.12b[) quotients out the order of red crossings in the succession of eras. Small billiard eras 
are defined in [lj] by the hit on the red wall, and a billiard map for this phenomenon has been pointed out to be much 
more 'intractable' than a CB-LKSKS map because the regions of the small billiard phase space are not defined by 
straight lines, and there exists no transformation able to map the curved domains into square boxes. This is consistent 
with the fact that the small billiard domain is not symmetric, in the sense that there exist no transformation able to 
identify its sides, as already pointed out. 
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The maps (|4.12l) avoids the problem of identifying sides of the small billiard by identifying the kind and the order of 
the 'crossings' of the symmetry lines, for which the definition of small billiard eras looks in principle lost. Anyhow, by 
using the map (|4.12j) . it is always possible to swich from the small billiard era description, where the hit on the red 
wall defines an era, to the quoticntcd big-billiard description, where BKL eras are described by the Poincare return 
map on the (identified) green wall(s). In other words, by the definition of a Poincare return map on the green walls, 
any information about the crossing of the symmetry lines of the billiard is disregarded. Anyhow, as shown in Table 
IIII1 this information can be recovered also in the BKL quotiented billiard. Also in the BKL-quoticnted billiard, the 
sub-regions are delimited by curved lines and are not box-like domains, and there exists no transfomation able to 
map these regions into boxes. 

From this point of view, one can infer also that, for the sake of a Poincare return map on the green walls, the kind 
and the order of crossings is unimportant. In fact, the r time evution of the scale factors defines the change of eras as 
the point where the non-oscillating scale factor inverts its slope. From the BKL map point of view, this point defines 
when two different sides have to be identified, while, from the CB-LKSKS map point of view, all sides are always 
identified, as described in Table llVl 

In the unquotiented dynamics, the change of era (i.e. the change of the corner where the oscillation takes place) is 
always accompanied with the crossing of a trajectory of a red symmetry line, and this crossing is not relevant for the 
point where the non-ascillating scale factor changes its slope, as one sees from Table IIIIl For this reason, an exact 
correspondence between the length of big billiard eras and small billiard eras is in principle lost, but can be recovered 
case by case, if needed, by gathering all the information together. 

k. More about the era-transition map on the UPHP By means of all these considerations, the presence of a 
(complex) reflection in the maps (|4.8[) and (|4.10[) can be further understood as follows. 

Following the example of a 3-epoch era of the ba type, followed by an ac epoch, two different subcases are possible. 
In particular, either the ac epoch is the first epoch of an ac era, i.e. according to the succession 

ba — > ab — > ba — > ac — > ca — > ... (4-13) 

or there is a 1-epoch era of the ac type followed by a cb epoch, i.e. according to the succession 

6a — S- ab — >• ba — >• ac — S- cb — S- ... (4-14) 

By construction, one can easily verify that, in the first case, i.e, according to the succession (|4.13[) . only one red 
symmetry line is crossed, i.e., in the small billiard version, the red wall is hit only once. Anyhow, the red symmetry 
line can be crossed either during the last ba epoch, or during the first ac epoch. 

On the contrary, in the case of the succession (|4.14[) , the red symmetry lines of the billiard can be crossed either one 
or three times. 

These considerations explain that 

1. there is a mismatch in the comparison of big-billiard eras and small-billiard eras, i.e. the red line can be crossed 
during different epochs; 

2. there is a mismatch also in the number of eras counted on the big billiard dynamics and on the small-billiard 
dynamics, i.e. the red symmetry line can be crossed even three times during the transition from different eras 
(recall that the crossing of a red symmetry line defines a small-billiard era) ; 

3. there is a (complex) reflection in the era-transition map on the UPHP, i.e. the number of crossings of the red 
symmetry lines corresponds to the number of complex transformations of the small-billiard dynamics; 

4. the role of the Kasner transformations in Table U can be interpreted as reflections with respect to the symmetry 
lines; the fact that half of them is complex and half of them is real is due to the further choice of re-orienting 
the type of eras with respect to the preferred ba type. 

I. Comparison of the maps on the UPHP According to all the steps followed in the above, a paradigm can 
be implemented, according to which it is possible to reconstruct all the maps starting from the knowledge of the 
desymmctrized fundamental domain, corresponding to the more general inhomogencous case, i.e. the geometrical 
symmctry-quoticnting. 

The paradigm can be summarized as follows 

1. take twice the desymmctrized domain, and obtain a symmetric domain with suitable boundary identification, 
corresponding to the fundamental domain of the group (|3.1[k 
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big billiard 


epoch map 


era map 


a — > a 
b^b 
c — > c 


a — > a 
b — > a 
c—>c 


a — > a 
b — ► a 

c — > a 



TABLE IV. Comparison of the unquotiented big- billiard map, the BKL epoch map, and the CB-LKSKS map (consisting of the 
composition of successive BKL era-transition maps), according to their different actions in quotienting the roles of the Kasner 
exponents, which is described by the different identifications of the boundaries of the billiard table. 

2. the unquotiented big billiard map is obtained by interpreting its domain as a suitable congruence subgroup of 
the desymmetrized domain without any boundary identification; 

3. the BKL epoch map for the big billiard is defined for a domain where only two sides are identified (in the ba 
case) by translation, and the other is identified with itself. The BKL epoch map for the big billiard is the 
generated by the itcration(s) of the identification; 

4. the BKL era-transition map for the big billiard is obtained by identifying all the sides of the billiard onto (a 
preferred) one, and is generated by (a suitable composition of) these identifications; 

5. the CB-LKSKS map for the big billiard is obtained by the iteration(s) of the BKL era-transition map; 

6. the CB-LKSKS map for the small billiard is obtained from that of the big billiard by inserting the reflection R3 
in the CB-LKSKS map where needed, and has the meaning of identifying the green sides of the big billiard by 
neglecting the information on the crossing of the symmetry lines; 

7. the maps for the reduced phase space (u~u + ) are obtained by evaluating the maps in the UPHP on the absolute 
v = 0. This maps act diagonally on each variable u~ and u + . 

On the other hand, the understanding of the BKL map in the (u~,u + ) plane allows one to understand the exact 
partition of the dynamics into epochs and eras by means of the 'boxes' B xy and F xy (given in Table II and Table 
III of [lH) and by means of the precise relations between the number of epochs in each era and the value of the 
variables lip, as given in [l4|, for which it is possible to work out the relation n& a — 1 = [u^ a ] used in (|2.5j) . The 
main features of these maps are summarized in Table IIV1 where the different symmetry quotienting mechanisms of 
the Kasner exponents are illustrated according to the different identifications of the sides of the billiard table. The 
three maps are defined on the same domain, given in Figure [TJ 



V. QUANTUM REGIME 

The quantum version of the billiard model is obtained by promoting to operators the momenta 7r a 's in the kinetic 
part of the Hamiltonian Hq =_JrG ?a *7r' 7r&, with the metric G a b defined in (|2.2|) . As a result [Hj], on the PUHP, the 
following eigenvalue problem [2(| is obtained 

- A LB *(z) =E9{z) (5.1) 

where Alb is the Laplace-Beltrami operator on the PUHP, and \l/ is the wavefunction solving the eigenvalue problem. 
From a cosmological point of view, the wavefunction \E' can be interpreted as the wavefunction of the universe, in the 
asymptotic regime when the universe approaches the cosmological singularity, and quantum effects have to be taken 
into account. 

The eigenfunctions VP can be analyzed according to their Fourier expansion for the periodic variable u [2fj| | . 
Maass wavefunctions f(z), whose Fourier expansion = $ s reads 

$ s (u,w;ci,c 2 ) = civ s +c 2 v 1 ~ s + ^c^y l/2 K, s _i{2-K \ fj, \ v)exp[2m^u], (5.2) 
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are the only (known) eigenfunctions for the Laplace-Bcltrami operator in (|5.1|) invariant under the modular group. 
The functions K, are the K-Bessel functions. The coefficient c M in (|5.2p can be imposed conditions for the well-behaved- 
ness of the wavefunction with respect to the polynomial growth of JC for z — > oo. The Maass wavefunction (|5.2j) then 
is defined by the energy level s and by the two constants c\ and C2. 
The eigenvalue problem (|5.ip for the Fourier-expanded Maass waveforms (|5.2[) is then 

- A LB $ s (u,u) = E s V a (u,v), (5.3) 

where the eigenvalues E s defined the energy spectrum of the quantum Hamiltonian, and explicitly read 

E s = s{s-1). (5.4) 

Maass wavefunctions solve the problem (|5.1|) . and obey the (invariant) transformation law f("fz) = f(z) V7 € 
SL(2,Z). 

Moreover, according to the features of the Fourier expansion for the u variable, they can be classified as u-odd and 
it-even, according to their parity under the transformation u — > —u. In particular, u-odd Maass wavefunctions contain 
only sin(27r/iu) in their Fourier expansion, while w-even one only cos(27r/iii). 

Finally, the behavior of Maass wavefunctions under reflections has been studied according to the parity under the 
reflection R\ (form which the behavior under reflections in general can be obtained). Accordingly, the Maass wave- 
functions can be classified according to their i?i-parity [22| . 

The quantum eigenvalue problem (|5.1[) will be analyzed from two different points of view: i) how to impose the proper 
boundary conditions according to the geometrical 'shape' of the billiard table; and ii) how to interpret the invariancc 
of the dynamics under the billiard maps. Analyzing the task in two different steps will be useful in further illustrating 
the differences between the geometrical symmetries of the billiard tables and the features of their dynamics. 
The following reasonings are in principle valid for both Dirichlet or Neumann boundary conditions. The choice of 
Dirichlet [l^j , [H| [33| |2l| or Neumann [23| [25[ boundary conditions for this particular problem and their comparison 
with respect to their quantum-mechanical interpretation has already been compared |28l | and discussed [22j in the 
literature. Anyhow, for the sake of compactness, only the Dirichlet problem will be mentioned; the way to treat the 
Neumann problem is analogous, and the comparison will be briefly discussed. 

m. Empty billiards Boundary conditions for the billiard table alone can be analyzed according to the 'shape' of 
the domains. 

For the 'empty' unquotiented big billiard, for Dirichlet boundary conditions, the eigenfunction $ must vanish on the 
three boundaries^, without any particular prescription on the dynamics. 
Automorphic forms f(z) that admit a Fourier expansion 

f( z ) = ^2 fn(v)exp[iiris] (5.5) 

n 

can be considered as eigenfunctions of the problem (|5.1[) . as their odd part would contain only the sin(n7ru) functions, 
while their even part only cos(n7rw), for Dirichlet or Neumann boundary conditions, respectively. 
We can therefore take for \P any automorphic form F, solution of (|5.1[) . vanishing, say, on u = 0, for the Dirichlet 
problem. In fact, without requesting any identification about the remaining two sides, it is straightforward to verify, 
from a geometrical point of view, that the line u = is mapped into the line u = — 1 (corresponding to the side b) 
by the transformation z — > z — 1 = T _1 (z), and that the line u = is mapped into the circle u(u — 1) + v 2 = 
(corresponding to the side c) by the transformation 2 — > —z/(z — 1) = STS(z), such that, given f(u = 0, v) = 0, W, 
then the proper Dirichlet boundary conditions follow, i.e. 

f{T-\u = 0, v)) = = f(STS(u = 0, v)) (5.6) 

These cigenfuntions can be further classified according to their behavior on the 'symmetry lines' of the big billiard. 

Boundary conditions on the small billiard table can also be imposed without any identification of the boundaries 
(which is indeed not possible), but according to simple considerations about the symmetry of the problem. 
In more detail, the appropriate wave functions for the 'empty' small billiard can be obtained by picking up from 
the u-odd wave functions of the 'empty' big billiard those that are antisymmetric (i.e. those that vanish) on the 
symmetry lines and that satisfy proper Dirichlet conditions. 



2 Although the eigenvalue problem 115.11 1 on the big-billiard table is well-defined, and although it is possible to treat it explicitly, in the 
literature other approaches are present, in which approximated or deformed domains are considered [3f| l'24f 
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As an example, we can choose, as in |26j . 

f(u, v) = C v Sm(w fJ,U)^/vK s (nUV), (5-7) 
v 

such that the Dirichlct boundary conditions for the empty big billiard hold, according to (|5.6I) . From these functions, 
one can further select those that vanish on the 'symmetry lines' of the big billiard, i.e. for [i = 2v, 

$ s (m, v) = Cfj, sin(2fiTTu)^/vK, s (2^TTv) = $ s (u, v, ; ci = c 2 = 0) : (5.8) 

the odd Maass wavefunctions (|5.2[) arc immediately obtained from (|5.2[) by choosing ci = c 2 = (for which cusp 
forms are defined), this choice being imposed by the request that the wavefunctions vanish on the boundaries. 

So far, one could be tempted to conclude that the wave functions of the small billiard consist only of a part 
(half) of those of the big billiard, when only boundary conditions are taken into account, according to the fact that 
the small billiard is a desymmetrized version of the big billiard, or, going the other way round, the big billiard consists 
of six (suitably glued) copies of the small billiard. However, the analysis of the dynamics will impose much stricter 
constraints, and will be able indeed to show their complete equivalence. 

n. The billiard ball After having analyzed the boundary conditions that can be imposed according to the geo- 
metrical 'shape' of the billiard tables, it is possible to investigate how the structures worked out in Section HVl impose 
further constraints. 

Taking into account the big billiard map 7~, the quotiented big-billiard maps and the small-billiard map t, one 
has to impose that the dynamics be invariant under this maps in any point, and not only on the boundaries. Further- 
more, the maps Ai on the UPHP impose further constraints on the realization of boundary conditions via the proper 
identifications of the sides. One is therefore forced to choose Maass wavefunctions (according to their properties of 
being invariant under the elements of the group SL(2,Z), to their it-parity and to their behavior under reflections) 
also for the big billiard both in the unquotiented case and in the quotiented ones. 
So far, distinctions have to be taken into account. 

The wavefunctions of the unquotiented big billiard have to satisfy suitable boundary conditions (without identification 
of the boundaries), and must have the appropriate behavior under the big-billiard map T (|4.2[) . 
The wavefunctions for the BKL epoch map must obey proper boundary conditions (with the identification of the 
boundaries b — > a and c — > c), and must be invariant under the transformation T" 1 . In fact, the BKL epoch map 
contains only an even number of reflections. Of course, this paradigm is valid only for a given era, i.e. for a finite 
number of epochs, and does not correspond to the description of the complete dynamics. 

The wavefunctions for the complete BKL map (i.e. when also the era-transition map is considered) and for the 
CB-LKSKS must obey suitable boundary conditions (with identification of the boundaries b — > a and c — > a), and 
must have the proper behavior under reflections. 

The wavefunctions for the small billiard can be obtained in principle, as explained in the discussion of boundary 
conditions alone, from those of the unquotiented case by 'picking up' those that vanish along the sides B and R (that 
can be interpreted as 'symmetry lines' for the big billiard). But, as already mentioned, the invariance under the maps 
forces one to consider only Maass wavefunctions for the unquotiented big billiard as well. As a result, because of the 
features of the Maass wavefunctions, the wavefunctions of the small billiard coincide exactly with those of the full 
unquotiented billiard. As the big-billiard map (|4.2p consists of a suitable composition of the small-billiard map (13.51) . 
the behavior under reflections will be the same. 

As a result, the wavefunctions of all the kinds of 4 = 3 + 1 cosmological billiards are the same, i.e. Maass waveforms. 

Maass wavefunctions can be classified as odd or even according to their behavior under it — > — u, i.e. according 
to their u-parity. The choice (|5.2I) corresponds then to it-odd Maass wavefunctions for Dirichlet boundary conditions. 
Furthermore, Maass wavefunctions can be classified also according to their behavior under reflections, i.e. R— 1-odd or 
i?i-even under the fundamental reflection R\ : z — > —z. The behavior under other reflections can be straightforward 
analyzed; as an example, the reflection i? 2 : z — > z" = — z + 1 can be analyzed as for z' = T~ 1 z = z — 1, such that 
z" = Ri(z'). i?i-odd Maass wavefunctions would therefore acquire a minus sign under the action of the billiard maps 
(which contain all an odd number of reflections). 

o. Comparison of the results As already pointed out, the invariance of the dynamics of the system under the 
action of the billiard map follows from the existence of a reduced symplectic form, whose integral over an appropriate 
region (of the (u~u + ) plane) accounts for the corresponding probability for a given sequence of epochs to happen. 
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From the quantum point of view, the probability is obtained by integrating over a suitable domain the squared 
modulus of the wave functions, for which a difference in sign does not produce any effect. More generally, physical 
observables are constructed from a suitable scalar product of wavefunctions. Accordingly, although the wavefunctions 
can have an even or odd behavior under reflections, this behavior is not relevant in the definition of probabilities [22 |, 
The same applies for u-parity. 

Similarly, as already pointed out in (22j | for the discussion of the choice of boundary conditions, Dirichlet and Neumann 
boundary conditions produce the same expression for the probability, as complex phases would cancel out. 
Nevertheless, it is worth remarking that a crucial difference between Dirichlet and Neumann boundary conditions 
arises. In the case of Neumann boundary conditions, i.e. for even Maass wavefunctions, also a continuum spectrum 
appears. 

Although the properties of the problem (|5.ip are already known, the comparison of the results for boundary conditions 
in the two different cases, i.e. the 'empty' billiard and the 'dynamical' billiard, is powerful in outlining the role of 
statistical billiard maps also in the quantum regime. In fact, as already pointed out, the billiard dynamics consists 
of a succession of free-flight Kasner trajectories joined at the reflections on the boundaries. In this evolution, the 
statistical maps encode all the information in the Poincare return map on (a preferred) side of the billiard. As a 
result, the continuous dynamics of the billiard is fully described by discrete maps. 

At the quantum level, it is then interesting to remark that the 'empty' boundary conditions account only for the 
geometrical shape of the billiard table, while the 'dynamical' boundary conditions (worked out of the statistical maps 
as the proper boundary identifications under the symmetry-quotienting mechanisms of the Kasner exponents) take 
into account the complete dynamics of the billiard (i.e. also the free-flight Kasner evolution) by imposing further 
constraints on the boundary conditions. 

Collecting all the ingredients together, it is possible to conclude that i) no complete symmetry-quotienting of 
the big billiard is able to eliminate reflections from the billiard dynamics. Anyhow, it is possible to identify in the 
BKL map two different symmetry-quotienting mechanisms, i.e. one that does not require the presence of reflections, 
and one that does. This difference implies a (formal) difference in the imposition of 'dynamical' boundary conditions, 
the difference being rather formal because of the invariance properties of Maass waveforms; ii) the choice of Dirichlet 
or Neumann boundary conditions does not affect the definition of probabilities, but the independence of probabilities 
from the behavior of Maass wavefunctions for i?i-parity and u-parity is somehow 'accidental', i.e. it is due to the 
explicit Fourier decomposition of the wavefunctions, and does not allow one in any manner to construct billiard maps 
which do not contain reflections (or which contain an even number of reflections); and Hi) considering symmetry 
(sub)groups generated by transformations which do not contain reflections can be anyhow useful to understand the 
the symmetries and the possible identifications of the sides for cosmological-billiard tables. 



Once the BKL map has been defined for the z variable on the UPHP, it is possible to determine subdomains on 
the big billiard table, that correspond to the first epoch of each BKL era, and to analyze the corresponding conserved 
quantities. 

The definition of the subdomains of the billiard table for BKL epochs and eras follows from transferring the informa- 
tion gained in the (u~,u + ) plane to the UPHP. 

In the meanwhile, the dynamics of each xy phenomenon covers the complete billiard table. The possibility to uniquely 
transfer the information to the UPHP is given by the BKL maps, as they allow one to choose a preferred (ba) kind 
of phenomena an to interpret them uniquely as fully representative of the unquotiented dynamics. Because of this, 
the definition of the subdomains of the billiard table by segmenting the u axis according to the BKL map is fully 
consistent. Nevertheless, these subdomains are overlapping and do not 'tile' the big billiard table uniformly. 



p. Recalling the features of the reduced phase space: the 'classical ' BKL probabilities In the reduced phase space 
(u~ , according to the paradigm implemented in [14j |. it is possible to define the probability P nicy for a n^-epoch 



era of the xy type to take place as the integral of the reduced form u in (|2.3I) over the corresponding box in F xy . 
According to the invariance properties of the reduced form w, and to the implementation of the BKL map, this 
probability equals the probability P n for a n-epoch BKL era to take place (expressed, for technical reasons, in Fb a )- 
The integration extrema can be defined independently for u + and u~ (as the integration domains are boxes), by 
selecting the range of u~ and u + corresponding to n. 

Furthermore, given a point u = (u~ , u + ) € -Bfc a , one can interpret the corresponding trajectory as the k-th epoch of a 
(n > fc)-epoch BKL era, whose probability to take place equals the BKL probability P n (by the invariance of the form 
uj under the BKL map). Moreover, the probability P ni .n 2 f° r a ni-epoch BKL era to take place and to be followed 
by a n2-cpoch BKL era composes by means of the definition of further sub-boxes in Fb a . 
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TABLE V. The boundaries of the subdomains, which correspond to a given BKL epoch, specified by iV. N — n denotes the 
first epoch of a BKL era, while N = n' < n denotes the n'-th epoch of a BKL era. 
Left panel. The boundaries of the boxes in the (u~ ,u + ) plane. 

Right panel. The corresponding values for the radius r and for the center uo of the geodesies parametrizing the BKL epoch 
in the UPHP. 



In particular, one can define the ranges of u and u + as a function of the features N of the epochs, i.e. 

u+(N) < u+(N) < u+(N), u-(N) < u-(N) < u M (N), (5.9) 

where the subscript m generically denotes the minimum value, while the subscript M generically denotes the maximum 
value. As a result, the (normalized) probability Pn is given by 

Pn = ^: du- du+— — . (5.10) 

The details of the ranges (|5.9[) for the three situations described in the above and their BKL probabilities are recalled 
in the Left Panel of Table [V] In the Right Panel of Table [V] for a given BKL era, the values of the radius r and of 
the center uq of the geodesies (|4.L|) parametrizing the corresponding BKL epoch arc listed (for future purposes). 
It is now useful to remark that the probability P n for an rt-epoch BKL era to take place corresponds to the area of the 
corresponding starting box on the (u~ ,u + ) plane, according to the measure l/(u + — u~) 2 . One needs now to identify 
the proper region of the UPHP ( i.e. the subdomain of the big-billiard table) which corresponds to the proper box of 
the (u~ ,u + ) plane. 

It is interesting to stress again that one of the most useful features of the reduced-phase space method is the possibility 
to define suitable regions and subrcgion in the (u~ , u + ) plane, which arc delimited by straight lines (boxes). According 
to this, the definition of the allowed subdomains for epochs and eras acquires a straightforward representation. On the 
contrary, on the UPHP, the definition of the corresponding subrcgions would appear technically more complicated. 
This way, the best approach is to transfer to the UPHP the information gained in the reduced phase space. By means 
of the BKL map on the UPHP, the definition of these subdomains will acquire a well-posed meaning. In fact, one 
should remember that each B xy region of the (u~,it + ) plane covers the complete billiard table. 

Transferring the information available in the reduced phase space, as listed Table [Vj one learns that these re- 
gion will be delimited by — 1 < u < and by the proper specifications for the trajectories of the billiard parametrized 
by geodesies G, as in (|4.ip . Following the notation of (|5.9|) . one then finds the integration extrema by plugging the 
results recalled in the Right Panel of Table fVl into (|4.1[) . In particular, the range v m (n) < v < v M (n) has to be 
worked out of (|4.1|) . 

One then finds that the domains D n correspond to the portion of the billiard table between two geodesies (|4.1[) Q M 
and Q m , defined as 

Qm = gM(u~ n ,uj f ), g m = g m {u~ hI ,ui n ). (5.11a) 

and the 'vertical' walls of the billiard, u = and u = — 1. The area a„ of each subdomain D n is evaluated according 
to the measure da on the UPHP, i.e. 

dudv 

da=^^, 5.12 

v 

such that the of the subdomain corresponding to the n-epoch BKL era is given by 

rv M (n) , 

du / % (5.13) 

-1 Jv m {n) 



As one verifies straightforward, the measure dudv / v 2 , i.e. the invariant measure on the UPHP, is invariant under the 
BKL epoch map and under the CB-LKSKS era transition map on the UPHP. As a result, the quantity a n is conserved 
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during the evolution of the billiard. 

One can also interpret (|5.13[) as the integral of the area form A(u, v), i.e. 

. . , du A dv . 
A{u,v) = 5 — . 5.14 

which is therefore a form invariant under the BKL map and the CB-LKSKS map. Its integral over suitable domains 
will provide conserved quantities. The invariant form (|5.14[) can also be considered to be obtained form the invari- 
ant ui form on the UPHP, when a suitable cross section of the Hamiltonian flow is taken into account, as stressed in [lij ]. 

The conserved quantities a n ( c ii c 2i s ) 

r ° . r M[n) , _ . . ,o dv 

-1 Jv m (n) v 



a n (cx,cr,s) = I du I | $ s (u,w;ci,c 2 ) | 2 (5.15) 



under the BKL epoch map and the CB-LKSKS map correspond to the regions of a BKL epoch in the UPHP, and 
are obtained by transferring on the UPHP the information gained in the restricted phase-space. They are obtained 
from the s-th Fourier mode of the decomposition <J> S of the wavefunction ^, and depend on the number n of the BKL 
map. It is interesting to remark that, from a classical point of view, one does not need to require that the coefficients 
ci and C2 In (|5.2j) and in (|5.15j) vanish. In fact, from the classical point of view, one only needs invariance under the 
BKL epoch map and the CB-LKSKS maps. The expression (|5.8[) is obtained, in the quantum regime, by imposing 
boundary conditions. 



VI. QUANTUM BKL MAPS 

One has now collected all the features of the quantum regime that allow one to express (at least formally) the 
'quantum BKL probability' for epochs and eras to take place. 

After the implementation of the quantum regime, the expression for the quantum probabilities, as the integral 
of the squared modulus of the Maass waveform over the suitable domain in the UPHP, follow directly. In fact, this 
has the direct physical interpretation of the quantum probability for the billiard ball to be in the corresponding sub- 
region. By means of the BKL map on the UPHP, and of the division of the billiard table into subdomains D n which 
correspond to an n-epoch BKL era, it is possible to obtain the quantum analogue of the classical BKL probabilities. 
Moreover, the invariance of the Maass wavefunctions under the BKL map will ensure the correct behavior of these 
probabilities during a BKL era. 



q. Quantum implementation of probabilities for BKL epochs and eras. As anticipated in the above, it is now 
possible to establish a formal definition for the quantum version of BKL probabilities for a n-epoch BKL era to take 
place, say P„ jS . 

As a result, the quantum BKL probability P„ will be given by the probability for the billiard ball to be found into the 
proper integration domain in the billiard table. In other words, one has to look for the probability that a geodesies 
Q (parametrizing a BKL epoch) be between the proper integration domain. 



r. Formal results Quantum BKL probabilities Pjv will then be given by the probability for the billiard ball to 
be found between the proper integration domain D n in the billiard table. They are expressed as the integral of the 
squared modulus of the Maass waveform (|5.8p over the proper integration domain D n , such that, in symbolic notation, 

Pn.s= / du / ^ | I 2 • (6.1) 

J-l Jv m (n) V* 

In (|6.ip . the wavefunction is $> s (u,v;ci = c 2 = 0) in (|5 .8[) rather than that of (|5.2|) . because the proper boundary 
conditions have to be recovered. 
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s. Physical interpretation of the quantum probabilities with respect to the BKL maps Given the formal result 
(|6.1[) , the physical interpretation of the quantum BKL probabilities can be understood. 

At the classical level, the reduced uj form is obtained by choosing ab initio a fixed energy shell, and then by 
considering the Poincare return map on the sides of the billiard. 

The BKL probabilities are defined as the integral of the reduced w form on the pertinent regions of the restricted phase 
space. Transferring the corresponding information from the restricted phase space to the UPHP, via the quoticnting 
of the roles of the Kasner exponents in the BKL map, one finds the corresponding subrcgions of the big billiard table 
by choosing the correct geodesies from their endpoints on the u-axis. 

On the other hand, the dependence of the quantum BKL probability on the energy level has therefore no straight- 
forward classical analogue. In fact, the eigenvalue problem (|5.1[) admits a complete spectrum for the energy levels 
E Sl and it is not possible to fix any value for the energy, and restrict the available phase space accordingly. On the 
contrary, it is however possible to define a quantum BKL probability for each energy level. As a result, the quantum 
BKL probabilities 'scan' the complete energy spectrum of E s and receive a contribution from each energy level s. This 
contribution is exactly the contribution of the generalized invariants (for a particular choice of the constants C\ and C2). 

Moreover, the quantum BKL probabilities (|6.ip inherit from the Fourier expansion of the Maass waveforms for 
the periodic variable u also the sum on the /i index of the u periods. As a result, the quantum BKL probabilities 
'project' the information about the length n of a BKL era along the directions /i of the it-periodic expansion. In other 
words, a sort of (generalized) Fourier expansions in the [i directions is obtained also for the values n, which defines 
the BKL length of an era. 

Given all these considerations, the quantum BKL probabilities are characterized, for each length n of the BKL 
eras, by 'sampling' all the energy levels of the eigenvalue problem (|5.ip via the corresponding generalized invariant, 
and by projecting them along the (i directions of the Fourier decomposition for the variable u. 



t. More about the UPHP Furthermore, given all the features of the quantum version of BKL probabilities, it is 
possible to comment on the generalized area invariants (|6.ip . 

In fact, as already mentioned, the generalized area forms do not owe their invariance to the choice of a fixed energy 
shell. The corresponding conserved quantities are defined for domains D n , which have been defined according to 
a fixed energy shell, to which uj owes its invariance in the reduced phase space. Accordingly, one can interpret 
the occurrence of all the s generalized forms as compensatory for all the possible classical values of the (classical) 
Hamiltonians, which have been integrated away to obtain w itself. 

On the other hand, from the classical point of view, the Maass functions (|5.2p are defined as invariant under the 
elements of the pertinent transformation group. In principle they can be investigated also without taking into account 
the eigenvalue problem; or, more in general, the eigenvalue problem can be investigated regardless to the meaning 
it acquires in quantum mechanics. But, in the quantum regime, the spectrum of the Laplace-Beltrami operator on 
the UPHP acquires the meaning of the succession of energy levels for a 'particle in a box' (i.e. the quantum billiard 
dynamics). This way, choosing real eigenvalues for the Schroedinger equation in the quantum regime implies the 
interpretation of the class of generalized area invariants on the UPHP as due to the real possible values for the 
(conserved) energy of the (classical) Hamiltonian, which are real. 



u. Open issues It is in order to remark that the definition of quantum BKL epochs and eras requires further 
investigation. 

On the one hand, the remormalizability of the integrals defining the quantum BKL properties on the UPHP, where 
suitable suregions of the billiard table have to be taken into account, seems to require numerical simulation. 
On the other hand, the role of the quantum BKL properties in the definition of the quantum nature of the cosmological 
singularity has not been clarified yet. In fact, in the pioneering paper by Misner [35(, the billiard table is examined on 
the unit circle rather than on the UPHP, and the quantum numbers that are hypothesized to describe the quantum 
nature of the cosmological singularity are obtained by considering a squared potential weel, raher than the exact 
traingular domain. Moreover, it is possible to demonstrate that the Misner variables correspond only asymptotically 
(towards the singularity) to the description on the UPHP. 

On the other hand, from the classical point of view, the BKL length of epochs and eras is crucial in the definition of 
the cahotic nature of the cosmological singurity, when an isotropic and inhomogenous universe is considered. In this 
respect, no previous interpretation of its quantum nature has been proposed. 
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Furthermore, the role of the energy shells defined by the Fourier modes p, of the decomposition of the wavefunction 
$ s has no classical analogue, and can be interpreted as the possibility to isolate only selected energy shells. In this 
respect, it is necessary to stress that the definition of BKL probabilities from the classical point of view is connected 
with the billiard description of the cosmological singularity on by fixing ab initio a particular energy shell in the 
Hamiltonian flow [3]; this suggests that the quantum definition of theese energy shells is not connected with the 
quantum description of billiards on the hyperbolic plane in general, but descends only from the invariance of the the 
proper integrals of the Maass waveform on the UPHP under the BKL map. Also in this direction, further investigation 
is needed to clarify the role of (a suitable selection of) energy shells in the description of the quantum nature of the 
universe. 

Nevertheless, the presence of these new structures in the wavefuntion on the UPHP has to be regerded to as the 
discovery of new structures in the quantum nature of the cosmological singularity. Their physical interpretation 
requires further developments also in view of its numerical computation. Indeed, the renormalizability properties of 
BKL probabilites in the classical framework depends crucially on the chioce of a fixed energy shell. In the quantum 
case, these properties has to be studied yet. Nonetheless, the case of a non-renormalizable probability set due to a 
non-standard measure is not to be considered as a flaw in the construction; in fact, also inn [14| . a non renormalizable 
measure on the unit circle has been defined. 



v. UPHP vs UH Several investigations in the literature face the problem of the quantum version of cosmological 
billiards on the UH. 

The full picture in which the Einstein field equations are mapped to is indeed that of the cahotic motion of a 
point particle inside the UH, which is coordinatized by the functions /3 M = p7 M . As already remarked, solving the 
Hamiltonian constraint for the gravitational action with respect to the variable A = In p is equivalent to projecting 
the motion inside the hypcrboloid onto the surface of the hypcrboloid itself (i.e. one degree of freedom is eliminated). 
After the elimination of the Hamiltonian constraint, the projection on the unit disk and on the UPHP are considered 
only for the sake of a better visualization and of aesier calculation. 

In [331 ]. the complete wavefunction (which contains also the variable p) is factorized, and, as a result, the complete 
wavefunction and all its p derivatives tend to zero as the cosmological singularity is approached, when Neumann 
boundary conditions are taken into account. More in general, the full wavefuntion is demonstrated to be in general 
complex and oscillating. 

Another advantage of using the projected model is the possibility to keep under control the initial data. In fact, all 
the initial conditions for which the continued-fraction decomposition of the initial values of u + and u~ for a given 
tra ject ory is limited (which is a set of non-zero measure in the space of possible data) are to be excluded [36| . 
In [37| , wavepackets for quantum cosmological billiards are proposed. The problem is envisaged also for the UPHP, 
where the wavepackets are constructed by means of a suitable Gaussian weighting function. The behavior of the 
model in the vicinity of the cosmological singularity is analyzed by computer simulation according to a very particular 
set of initial data. 



VII. CONCLUDING REMARKS 

In this paper, we have analyzed he group structure underlying cosmological billiards in 4 = 1 + 3 dimensions. In 
particular, we have taken into account the unquoticnted big billiard, the ^-symmetry quotiented billiard, and the 
small billiard. The BKL map and the CB-LKSK map consist in quoticnting out of the dynamics the (6 possible) 
permutations of the Kasner exponents, while the small billiard is obtained in the inhomogeneous case, and consists 
namely of (1/6) of the big billiard. 

The implementation of the quantum regime has been carried out on the basis of both boundary conditions and 
properties of the dynamics. As a result, a difference in the imposition of the boundary conditions appears, but this 
difference is eliminated as soon as the dynamics is taken into account. More in detail, although a more general 
wavefunction can be chosen in the case of the unquotiented big billiard as far as boundary conditions are concerned, 
the Maass wavef unctions remains the only possibility when the dynamics has to be taken into account. 
An investigation of conserved quantities under the BKL map and the CB-LKSKS map in the UPHP at the classical 
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level allows one to introduce a quantum analogue for the BKL probabilities. The quantum version of the BKL prob- 
abilities for the length n of an era has been established as the integral of the squared modulus of the eigenfunction 
of quantum the Hamiltonian over a suitable subdomain of the big-billiard table, and can be compared with the 
occurence, at a classical elvel, of generalized area forms on the UPHP. As a result, the BKL probabilities experience 
the energy levels s of the spectrum of the quantum Hamiltonian and project each contribution along the directions fi 
of the Fourier series of the Maass wavefunctions for the periodic variable u. 
The work has been organized as follows. 

In Section HI the main issues about cosmological billiards have been introduced. 

In Section [Til 4 = 3 + 1 cosmological billiards have been obtained, and the problem of suitably quoticnting the 
dynamics has been focused on. 

In Section IIII1 the main features of some groups have been recalled, as they prove to be a very powerful tool in 

describing the dynamics of the complete big billiard, as well as those of its symmetry-quotiented versions. 

In Section IIV1 the properties of the BKL epoch map, the BKL era-transition map and the CB-LKSKS map are 

compared. Although for these maps the properties of the dynamics are encoded in a description that destroys the 

typical behavior of the billiard, consisting of bounces (reflections) against the sides of the billiard table, a difference 

is found between the BKL epoch map and the CB-LKSKS map, interpreted as a composition of BKL era-transition 

maps. This is due to the different approaches in quoticnting the roles of the Kasner exponents. 

In Section [Vj particular attention has been devoted to the definition of invariant quantites on the UPHP. 

In Section [VTl the quantum implementation of the model has been addressed. More precisely, prescriptions about the 

the quantization are obtained considering both boundary conditions and maps. 

In Section IVTI a particular case of observables has been considered: the quantum versions of BKL probabilities for the 
length n of a BKL era. The quantum BKL probabilities for each n have been shown to depend also on the particular 
energy level s of the spectrum of the quantum Hamiltonian, and to project this information along the directions p, 
of the Fourier expansion of the Maass waveforms for the periodic variable u. The main results are summarized in 
Section [VlTJ 
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FIG. 2. The small billiard table on the Poincare upper half plane is sketched (thick solid black line). 

The 'symmetry lines' of the big billiard (corresponding to the other subdominant symmetry walls for the inhomogeneous case) 
are also displayed; in particular, the blue (dashed) lines are the part of the symmetry walls that bisects a given corner, while 
the red (dotted) lines are the part of the symmetry walls perpendicular to the side opposite to the given corner. The green 
(solid thin) lines are the sides of the big billiard (which correspond to the other curvature walls). 

The colors (dotted lines, dashed lines and thin solid lines) of the symmetry lines of the billiard describe how to 'glue' the 6 
equivalent copies of the small-billiard table and obtain the big-billiard table in Figure [1] In more detail, the G(reen) side has 
to coincide with the green (solid thin) line, the B(lue) side has to coincide with the blue (dashed) line, and the _R(ed) side has 
to coincide with the red (dotted) line. 
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FIG. 3. The difTerent regions of the (u~vA) plane considered in the different billiard maps. The (thick solid) black lines delimit 
the boxes B xy for the epoch hopscotch of the big billiard. The (dotted) blue lines delimit the starting boxes F xy of the era 
hopscotch for the big billiard, their details are given in The (dashed) purple lines delimit the domain where the dynamics 
of the small billiard takes place. 
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FIG. 4. Integrated description of the UPHP, which shows the big billiard table, the small billiard table, the symmetry lines of 
the big billiard wich define the small billiard, and an oriented geodesies that defines a trajectoty in the billiard. 
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FIG. 5. The subregions of the phase space of the Fba box of the billiard hopsotch court, where the BKL map acquires different 
features with respect to its content of Weyl reflection. The functions u a (green), up (yellow) and « 7 (blue) are plotted. The 
Fba box is delimited by the orange dotted line. The value u~ — — $ (red dashed line) divides the F oa starting box, as explained 
in Table ITlTI and in the text. 



